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EXPERIMENTAL STUDY BUBBLES MOVING LIQUIDS 


SYNOPSIS 


connection with other investigations the David Taylor Model Basin 
fundamental study the motion bubbles was undertaken. initial 
step, experiments were conducted determine the drag and shape single 
air bubbles rising freely various liquids. 

The results the experiments show that complete description the 
motion air bubbles use three dimensionless parameters containing 
the usual physical properties the liquid (viscosity, surface tension, density) 
not possible. Three types bubble shapes were observed each liquid, 
namely spherical, ellipsoidal, and spherical cap. For specific liquid the 
shape the bubble was function its volume. 

For tiny spherical bubbles the drag coefficients coincide with those 
corresponding rigid spheres. With increase bubble size, decrease the 
drag compared that rigid spheres occurs some liquids. Thus, the 
drag curves the spherical bubbles rising various liquids fall between 
two limiting curves, namely the drag curve rigid and fluid spheres, respec- 
tively. was not possible determine criterion for the transition the 
bubbles from rigid fluid spheres. Ellipsoidal bubbles occur different 
ranges Reynolds numbers for the various liquids. The drag coefficients 
spherical cap bubbles are independent bubble size and have constant val- 
2.6. 

For bubbles (equivalent radius 0.03 0.30 cm) rising tap water 
water containing certain surface-active substances, experiments show 
increase drag compared bubbles pure water. 

Results tests determine the effect the container walls the ve- 
locity rise are also presented and description the experimental appa- 
ratus given. 


INTRODUCTION 


detailed knowledge the mechanism bubble motion interest 
the many fields engineering which deal with the prediction performance 
equipment and processes involving the motion gaseous bodies. 
general, the equation motion such bodies given by: 


Drag Pressure Force Weight (Mass Added Mass) Acceleration 


This equation can solved the total mass the bubble and all forces 
acting are known. 

The drag bubble is, general, complicated function its geometry, 
its velocity, and the physical properties the medium. The shape which the 


Physicist, David Taylor Model Basin, U.S. Dept. Navy, Wash., D.C. 
Mathematician, David Taylor Model Basin, U.S. Dept. Navy, Wash., D.C. 
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bubble assumes turn complex function the hydrodynamic, viscous, 
and interfacial forces exerted. addition, such effects container walls 
and bottom, free surface the liquid and the effect adjacent bubbles may 
also have pronounced influence the drag bubble. Unfortunately, 
there only meager information available the drag gas bubbles and 
practically none their added masses. 

initial obtain information the drag gas bubbles, study 
single bubbles rising freely their terminal velocity various liquids 
was undertaken. For the case motion such constant velocity and above 
equation reduces to: 


Drag Buoyant Force Weight 


since all forces acting the bubble are equilibrium. Specifically, the 
experiments described this paper consisted the determination the 
terminal velocity, shape, and path single air bubbles rising freely vari- 
ous liquids function bubble size. The possible effect the walls 
the container the velocity rise was also investigated. Future studies 
the effect different pressure gradients the drag and studies apparent 
masses gas bubbles are contemplated. 

Theoretical solutions for the drag bodies, applicable gas bubbles, 
have been obtained for spheres moving slowly infinite medium. 


For rigid spheres: Stokes (1) 


where the drag the sphere, the coefficient dynamic viscosity 
the liquid medium, the radius the sphere, the terminal velocity the 
sphere, and the coefficient surface For the fluid spheres (in 
contrast the rigid case) non-zero tangential velocities the sphere inter- 
face are assumed. Hadamard and Rybczynski assumed equality the tan- 
gential velocities and stresses the bubble-liquid interface, while Boussinesq 
included changes the interfacial stresses due ‘dynamic’ surface tension. 

the last the above equations, the factor approaches for 
approaching zero very large; approaches 2/3 for i.e. for 
large for approaching zero. Hence, for very small bubbles, Boussinesq’s 
solution approaches Stokes’ law limit, while the other limit the 
Hadamard-Rybczynski solution. 

dimensional analysis the motion rising bubbles yields when con- 
sidering eight physical variables3; 


(Cy R, W, H/y' , P/p') = 0 


These physical variables are: the velocity the bubble (U), the 
acceleration due gravity (g), the density the liquid medium 
the density the bubble the length parameter the bubble 
the coefficient dynamic viscosity the liquid and the gas 
and the surface tension (¢). 
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the Weber number the Reynolds number inside the bub- 


sity and viscosity the gas are considered negligible, the number dimen- 
sionless groups are reduced three: 


etc. For bubbles, the most convenient length parameter the equivalent 


radius defined The drag coefficient then becomes 


ficients are, respectively 
sphere) 


For the Stokes and Hadamard-Rybczynski solutions the drag coef- 


For the special case which only the velocity, the acceleration gravity, 
the density the liquid, and the equivalent radius are considered, one dimen- 
sionless parameter obtained, namely constant. This solution will 


shown apply the region very large (spherical cap) bubbles. 

Previous experimental work the rising motion bubbles has been ex- 
tensive. The investigations have, however, largely been restricted bubbles 
rising water. detailed bibliography these investigations may found 
Reference (4). 


Experimentation 


The experimental study consisted measuring the terminal velocity 
individual bubbles various sizes rising eight liquids. also included 
the determination the effect the walls the container the bubble 
velocity. Details the experimental apparatus and procedure, the genera- 
tion the bubbles, and the test liquids are given the following paragraphs. 

Test Tanks and Liquids. The tests were performed three transparent 
wall tanks; the large one was cross section and ft. height, the 
medium one ft. cross section and 1/2 ft. height, while the small one 
was inches cross section and ft. height. addition, tests were also 
performed insert inches cross section and inch height placed 
the center the medium tank. 

Since the large tank was sufficiently large dimensions, appreciable 
wall effect was expected.* The medium tank was chosen dimensions large 


few the previous TMB tests (5) were repeated the large tank 
observe any change results. These previous tests were conducted ina 
tank 1/2 ft. cross section and ft. height with ft. depth fil- 
tered water room temperature, using one end the tank for the tests. 
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enough reduce wall effects, yet small enough allow use variety 
liquids. The small tank and insert provided additional tank size. was 
intended, wall effect existed, extrapolate the results obtained the fi- 
nite containers the case infinite medium. 

The eight test liquids were water (at three different 
methyl alcohol, water containing 0.42% (by volume) Glim®, turpentine, mineral 
oil, two corn mixtures. Turpentine was selected one 
the test liquids because has, room temper?ture, the same viscosity 
cold water. One the corn syrup mixtures had approximately the same vis- 
cosity the mineral oil. The viscosity the liquids was measured means 
ordinary and modified Ostwald The accuracy measure- 
ment viscosity was 1/2% and 1/2%, respectively. The surface tension 
was determined the capillary-rise method (accuracy measurement: 3%) 
and the specific gravity the liquids was obtained means hydrometers 
(accuracy measurement: 1%). These physical properties were measured 
following the completion each test. They are summarized Table to- 
gether with those liquids used several other investigators. (6, 8). 

Varsol and water (room temperature and hot), tests were conducted 
all three tanks; for cold water and mineral oil, the medium tank and insert; 
and for all other liquids, the small tank only. 

Test Procedure. Small bubbles were generated means hypodermic 
needles and glass nozzles various sizes. The larger bubbles were obtained 
use dumping cup, which was inverted release the air bubble. The 
nozzles and needles were connected brass tube which was fastened 
sliding mechanism (Figure 1). This sliding mechanism allowed the tips 
the various nozzles placed the identical position, eliminating the need 
for refocusing the camera after each change nozzle. The air was supplied 
from compressed air bottle. needle valve regulated the air flow that 
bubbles were released the interval desired. 

The bubble size was determined ‘weighing’ sufficient number bub- 
bles the inverted funnel (Figure means analytical balance. Since 
the density air negligible comparison that liquid, the difference 
balance reading equals the buoyancy the bubbles (i.e. equals the vol- 
ume the bubbles times the density the liquid). The change balance 
reading was always least 0.2 gm, resulting accuracy measurement 
volume 1%. The volume the individual bubble was obtained divid- 
ing the total volume the number bubbles collected the funnel. com- 
parison photographs different bubbles showed that the bubble size did not 
vary the frequency bubble generation remained constant. Large bubbles 
from the dumping cup were weighed individually. The volume the bubble 
was adjusted for the change pressure due difference depth between 
the level which the rate rise determined and the level the inverted 
funnel. This was done use the general gas law constant temperature 
taking into account the partial pressure the saturated vapor test 
perature.? Tiny spherical bubbles could not generated frequency 


trade name (Standard Oil Company) for mineral spirits (heavy naphtha), 
petroleum derivative. 

Glim, surface active agent, the trade name (B. Babbitt, Inc) 
non-ionic, liquid detergent, condensation product ethylene oxide and 

Details this correction may found References (4) (5). 
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allow sufficient number collected the funnel, hence their size was 
determined from the photographic record. correction for change depth 
then needed. 

avoid any changes the volume the bubble due air interchange 
with liquid, the latter was saturated with air prior actual testing. This 
was accomplished stirring the liquid and blowing air through it. The 
liquids for all but the cold and water tests were room temperature, 
which varied little throughout the day. Water was cooled circulation 
through water cooler; was heated immersion heaters obtained di- 
rectly from the hot water faucet. Both filtered and tap water were used 
the tests. 

Uniformity liquid temperature was achieved means mechanical 
stirring before each test. Frequent checks temperature various loca- 
tions inside the tank were made means immersion thermometers. 
the process stirring, small bubbles appeared the liquid. The irregu- 
larity motion these small bubbles, which were still present after com- 
pletion the stirring, served indication the presence residual 
turbulence the liquid. sufficient time, the motion the small bubbles 
always became regular and hence indicated that the residual turbulence, 
still present, was not large enough affect the motion the bubbles. The 
actual test was not begun until all the small bubbles reached the surface 
the liquid. 

The rate bubble flow was then regulated the needle valve 
release bubbles with minimum spacing inches. This reduced the ef- 
fect the wake created the passage bubble the motion bubble 
following. The same precaution was observed for the larger bubbles that 
were formed dumping. additional precaytion was rotate the dumping 
cup with steady speed order avoid splitting the bubble the formation 
satellites upon release. The slow passage the air through the brass tube 
inside the tank allowed the air reach the temperature the liquid. Contact 
the air the nozzle tip inside the dumping cup with the liquid allowed 
saturation the air with liquid vapor, that the air bubble can, each 
instance, assumed saturated with the vapor the liquid which 
rises. 

Motion Pictures and Their Evaluation. The velocity, path, and shape 
the bubbles were obtained from motion pictures made with Mitchell 35mm 
‘camera special lens attachment permit close-ups. Film speeds 
frames per second and back lighting from white reflector were 
used. For the first few tests, the film speed was obtained photographing 
rotating clock dial; subsequently neon timing light (located inside the 
camera) with 60-cycle voltage source was utilized. From the marks the 
timing light the film, the film speed was determined. The field the 
camera ranged from 1.4 1.8 1.75 2.3 inches depending upon the refrac- 
tive index and horizontal depth the liquid. transparent scale photographed 
the plane the bubble provided the distance scale factor for the evaluation 
displacement and size. The camera lens was placed approximately the 
midpoint between liquid level and bottom tank for all tests. The camera 
location was each instance located sufficiently above the nozzle tip that 
the bubbles reached their terminal velocity before passing front the 
camera. 

Changes bubble volume due differences liquid depth were mini- 
mized making velocity measurements over very short vertical displace- 
ment (less than 1/2 inches). The rate rise bubbles was determined 
measuring the displacement bubble from reference point successive 
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frames the film means Bausch and Lomb contour measuring pro- 
jector using magnification twenty-five. These displacements were then 
plotted against the frame number. straight-line plot indicated that the ve- 
locity the bubble remained constant during the time passed the field 
the camera. From the slope the line, the frame speed, and scale factor, 
the velocity the bubble was computed. 


Results 


Terminal Velocity Bubbles. The results tests determine the ve- 
locity rise gas bubbles various liquids are most conveniently presented 
function the equivalent radius the bubble, defined the radius 
sphere having the same volume the bubble. Typical results are shown 
Figures and for filtered water (including data from other investigators) 
and for tap (unfiltered) water. summary all TMB velocity curves (ex- 
cept those for tap water) and those obtained Bryn given Figure 
compilation the properties the liquids was given Table general, 
the results seen from Figure indicate that for small (spherical) bubbles 
given volume the viscosity the liquid the most important property 
determining the rate rise. Very large bubbles (spherical caps) rise inde- 
pendently the properties the liquid. 

The results tests determine the effect the container walls the 
velocity bubble, conducted tanks different sizes water, Varsol, 
and mineral oil, showed within the limits experimental accuracy the ab- 
sence any wall effect for the range bubble sizes tested. Figures and 
are typical these results. Subsequent tests the other liquids were made 
the small tank only and its results may applied the case infinite 
medium. 

systematic investigation vertical proximity effect was made, i.e. the 
effect the wake created the passage bubble the motion bubble 
rising distance below. Such effects were avoided the experiments 
sufficiently spacing the bubbles. However, the results few special obser- 
vations indicate that proximity effects may appreciable. For 
tests mineral oil show increase and 39% for bubbles equivalent 
radius 0.17 cm, rising 7.7 and 3.2 apart respectively. The presence 
the wake the liquid thus results higher velocities rise the bubble. 

Figure compares the rate rise air and oxygen bubbles distilled 
water. significant change the rate rise the bubbles with change 
the gas composing indicated. is, therefore, reasonable assume that 
the velocity curves for bubbles composed various other gases would 
identical those for air bubbles. 

Nondimensional Evaluation. The results the TMB bubble tests and 
those Arnold, Bond and Newton, and Bryn were plotted terms three 
dimensionless products. Figure presents the results terms drag coef- 
ficient, Reynolds number, and the parameter ‘M’, while Figure presents 
them terms the drag coefficient, the Weber number, and The 
curve for filtered distilled water temperature 19°C was drawn 
through points obtained from the experiments Bryn and the TMB tests. 
Examination Figure shows systematic arrangement the curves 


The results for tap water and for water containing Glim are not shown. 


They will discussed subsequent sections. 
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with change the parameter which constant for specific liquid. 
can, therefore, concluded that neither the nondimensional sets presented 
nor any other complete set using the same six variables (namely velocity, 
acceleration gravity, density and viscosity the liquid, surface tension, 
and equivalent radius) sufficient for complete description bubble motion. 

The question now arises whether correlation bubble data could ob- 
tained using two additional dimensionless parameters, for example the 
liquid air viscosity and density ratios the Reynolds number inside the 
bubble and the density ratio etc. The results the experiments conducted 
not permit conclusions regarding the importance these parameters. 
short discussion the significance the internal Reynolds number will 
given subsequent sections. 

Three types bubble shapes were observed the experiments. Very 
small bubbles are spherical. Larger bubbles are flattened, i.e. ellipsoidal 
shape, whereas very large bubbles assume spherical cap shape. course, 
the volumes which these transitions occur vary with the liquids. Photo- 
graphs typical shapes are shown Figure should noted here that 
some the shapes shown these photographs are instantaneous shapes, 
since the shape large bubbles does not remain constant during the ascent, 
except for those bubbles rising highly viscous medium (e.g. mineral oil 
and corn cyrup). 

Table shows the approximate extent the three shape regions afunc- 
tion Reynolds number for the various liquids tested TMB. From this 
table and Figure (which also includes the drag curve for rigid spheres (11) 
and the lines for the Stokes and the Hadamard-Rybczynski law) the following 
can observed: The drag curves spherical bubbles the various liquids 
fall between two limiting curves. upper limit, the drag curve rigid 
spheres obtained, while the lower limit the drag curve for fluid spheres. 
With decreasing Reynolds number the rigid sphere curve connects with the 
straight line Stokes’ Law, while the fluid sphere curve connects with the 
line Hadamard-Rybczynski’s Law. The curve for the fluid spheres was 
obtained drawing the lower envelope the experimental curves; its ac- 
curacy can confirmed additional tests other liquids extension 
the theoretical solution into regions beyond that very slow flow. 

will also noted from Figure that the curve for mineral oil, for 
example, follows the straight line Hadamard-Rybczynski’s law over cer- 
tain region Reynolds numbers. This indicates that the boundary conditions 
assumed the analytical solutions for fluid spheres are actually fulfilled and 
that circulation exists inside the bubble. The experimentai curves also 
indicate interesting aspect the phenomenon bubble motion, namely 
that with decreasing Reynolds number, the drag coefficient the bubbles 
becomes equal the drag rigid spheres. The transition may occur 
Reynolds number about (as for filtered and distilled water) may not 
take place until very low Reynolds numbers are reached, i.e. well within the 
region slow flow (as for olive oil (7) very viscous syrup (6) Thus, 
from the experimental data available appears certain that tiny air bubbles 
rising any liquid follow Stokes’ Law. 

For bubbles behaving like rigid bodies, thus indicating absence motion 
inside the bubble, the internal Reynolds number (although non-vanishing) 
significance describing the rising motion the bubbles. Likewise, 
the internal Reynolds number cannot used predict the transition point 
which the drag the bubbles becomes less than that corresponding rig- 
spheres. Beyond this transition point, the internal Reynolds number might 
importance the description the motion the bubbles, 
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Surface tension tends make the surface area the bubble small 
possible. For given volume, the configuration minimum surface area 
This effect surface tension would most pronounced for bub- 
bles small radii, For larger bubble sizes, the surface forces become 
smaller compar.son the viscous and hydrodynamic forces and flattening 
the bubble occurs. This flattening approximately oblate spheroid re- 
sults higher drag compared sphere the same volume. From 
Figure and Table will observed that the region ellipsoidal bubbles 
for the various liquids occurs different ranges Reynolds number; that 
for liquids low ‘M’ number minimum the drag curves reached 
Reynolds numbers the order 250; and that these minima occur near the 
transition from spherical ellipsoidal shape. Such minima are not obtained 
for liquids high ‘M’ number. The drag coefficients bubbles such liq- 
uids decrease until constant value for the drag coefficient (spherical caps) 
attained. the ellipsoidal region the curves are arranged according 
the magnitude the ‘M’ number, indicating that the liquid properties con- 
tained this parameter, namely surface tension, viscosity, and density, are 
primary importance the motion these bubbles. Thus, this region 
the Reynolds number inside the bubble importance the description 
bubble rise, since correlation was obtained terms drag coefficient, 
Reynolds number, and ‘M’ number. 

When the viscous and surface tension forces become small relative the 
hydrodynamic forces, the shape assumed the bubbles that the so- 
called spherical caps. Typical shapes these bubbles are shown Figure 
The upper surface essentially spherical, while the lower surface varies 
from highly irregular one for liquids low viscosity smooth surface 
for very viscous liquids. The configuration the upper surface results al- 
most exclusively from the hydrodynamic The results the pres- 
ent tests number liquids confirm the constant value 2.6 for the drag 
coefficient obtained previous investigators (5, 12). The velocity ofspheri- 
cal cap bubbles given size rising any liquid can determined from the 
value the drag coefficient directly from the velocity curve (Figure 4). 
For 2.6, obtain for the rate rise the spherical 


caps all liquids 


Path Rising Bubbles. Figure shows typical paths and corresponding 
shapes bubbles four liquids. Three types motion the bubbles were 
observed the experiments: (1) rectilinear motion (2) motion helical 
path and (3) rectilinear motion with rocking. The motion spherical bub- 
bles either rectilinear helical. For ellipsoidal and spherical cap bub- 
bles all three types motion can occur. appears that the type motion 
may predicted from the value the Reynolds number which the motion 
takes place. Below Reynolds numbers about 300 the motion rectilinear. 
With increase Reynolds number spiraling begins and increases ampli- 
tude and frequency until maximum reached. Reynolds numbers 
about 3000 the spiraling disappears and only rectilinear motion with rocking 
obtained. For the bubbles rising mineral oil and the corn syrup mix- 
tures, only rectilinear motion without rocking was observed. The maximum 
Reynolds number reached during those tests was 150. 


See e.g. Reference (12). 
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The helical path the bubbles assumed either clockwise counter- 
clockwise direction, depending upon conditions generation. The velocity 
rise these bubbles not the sense which the bubble re- 
volves. The major axis ellipsoidal bubbles always perpendicular the 
direction motion. 

The oscillatory motion bubbles probably caused the periodic shed- 
ding vortices behind the bubble. Such vortex shedding has been observed 
experimentally for rigid spheres the same magnitudes Reynolds num- 
bers for bubbles (13). 

Bubbles Tap Water and Water Containing Surface-Active Substances. 
The drag coefficients for air bubbles rising tap water two different tem- 
peratures are given Figure 10. Gorodetskaya’s (14) results tap water 
room temperature are also included. comparison, the drag curves for 
bubbles filtered and distilled water room temperature and for rigid 
spheres are also shown. the region spherical and spherical cap bubbles, 
the drag curves the two temperatures coincide. The value the minimum 
drag coefficient is, however, greater than that the corresponding filtered 
water. general, for Reynolds numbers about 300, the drag curves 
bubbles tap water follow closely the curve rigid spheres.10 Thus, the 
results the experiments show that for bubbles (ranging equivalent radius 
from 0.035 0.25 cm) important whether the motion occurs filtered 
(distilled) tap water. 

view the fact that merely filtering the water was sufficient produce 
change the drag the bubbles, indicated that the presence minute 
particles causes this change. Most these particles, which are known 
exist ordinary tap water, can removed filtering. Specifically, such 
particles are present, high concentration them would found the sur- 
face the bubble. The particles the surface would travel with the bubble, 
hence imparting, effect, rigid surface the bubble. With increasing bub- 
ble velocities, the shear forces become large comparison the forces 
holding the particles the surface and hence certain critical velocity 
high concentration particles the surface can exist. 

The effect surface-active substances the rate rise air bubbles 
has previously been investigated Gorodetskaya (14), who added small con- 
centrations various alcohols water and concluded that, beyond certain 
critical concentration the surface-active substance, the rate rise the 
air bubbles not affected. Stuke (9) ran experiments with oxygen bubbles 
rising water containing small concentrations caproic acid. The concen- 
trations the alcohols and the caproic acid were relatively small, hence the 
decrease the surface tension was only about dyne/cm. The authors con- 
ducted tests water containing Glim, liquid detergent. the present 
study the concentration Glim (0.42% volume) was high enough decrease 
the surface tension dynes/cm. measurable change the viscosity 
and density the test liquid due the presence Glim was noted (see 
Table 1). This was also true for the alcohol and caproic acid solutions. Re- 
sults from these experiments are presented terms the drag coefficient 
and Reynolds number Figure 11. 

The drag curve for bubbles the Glim solution, well the experi- 
mental data from the other investigators water containing least the 


10. should noted here that the physical properties tap water did not 
differ from those filtered water. 
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critical concentration the surface-active substance, follows the drag curve 
rigid spheres Reynolds number about 200.11 the region Rey- 
nolds numbers 200, the drag curve for bubbles rising pure liq- 
uid having number very close that the Glim follows the 
drag curve fluid spheres. Thus, the motion bubbles water containing 
surface-active materials cannot compared with that bubbles pure 
liquids the basis drag coefficient, Reynolds number, and number, 
even the region ellipsoidal bubbles. Although the range bubble sizes 
did not extend fully into the region spherical cap bubbies, quite cer- 
tain that the presence surface-active substances will not alter their rate 
rise, which was shown independent all physical properties the 
liquid. 

The difference behavior the bubbles must sought the behavior 
the surface. high concentration molecules the surface-active sub- 
stance will found the surface the bubble. for the case the tap 
water, these molecules would travel with the bubble and impart, effect, 
‘rigid’ surface the bubble, that say impose the condition zero ve- 
locity the boundary. 

Thus, the results the tests given Figure show that these surface- 
active substances increase the drag the bubbles (in the region bubbles 
having equivalent radii 0.03 0.30 cm); beyond the critical concentrations, 
any increase concentration has relatively little influence ti.e drag the 
bubbles. 

Bubbles Rigid Bodies. previous sections was shown bubbles 
rising pure liquids behave essentially like fluid bodies over 
range bubble size, but that below certain critical size (the size 
different for various liquids) the bubbles behave like rigid bodies, that 
say the drag the bubbles equals that corresponding rigid bodies. 

possible explanation the anomaly behavior the gas bubbles 
given Boussinesq’s dynamic surface tension. assumed that inter- 
faces motion dynamic surface tension exists. Its magnitude given 
the sum the usual (static) surface tension and the dynamic increment. 
This increment includes constant proportionality (surface viscosity) 
which function the two fluids composing the interface. pointed out 
the ‘Introduction’, for small bubbles the effect the dynamic increment 
increases the drag the value corresponding rigid bodies. With increase 
bubble size this effect becomes negligible and the drag the bubble equals 
that fluid body. Therefore, for bubbles the transition region from ‘rigid’ 
fluid bodies would different for the various liquids. There is, however, 
experimental evidence that dynamic surface tension, postulated 
exists. 

From ‘hydrodynamic’ point view the reason for the transition the 
bubbles ‘rigid’ bodies not clear. The mere surface tension 
pressure drop the boundary conditions does not alter the 
analytical solution for fluid spheres. Hence, appears that the presence 
surface tension should have effect the motion the bubble, except 


11. For concentrations below the critical, the drag curve lies between that 
pure water and the curve shown Figure (9). 

12. The ‘M’ number Glim was 2.78 that Varsol, for example, 
was 4.3 

13. Mixtures, such the 13% ethyl alcohol-water mixture are included 
this category. 
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maintaining the spherical shape. Thus the anomalous behavior the bubbles 
cannot explained terms ‘hydrodynamics’, but must sought terms 
surface phenomenon. could shown analytically (if only the re- 
gion slow flow) that equality drag corresponding rigid bodies and bub- 
bles also implies equality boundary conditions the surface, then, 
the case rigid bodies, the velocity the entire surface the bubble must 
vanish. The surface must then able hold molecules the pure liquid, 
just the case tap water and surface-active substances the surface 
attracts and holds high concentration particles molecules the 
surface-active substance. The molecules the surface would travel with 
the bubble and hence would, effect, give the same boundary conditions 

rigid surface. However, the shear forces become larger comparison 
the forces holding the molecules the surface, ‘rigidity’ the surface 
cannot maintained; circulation inside the bubble ensues and the drag 

the bubble becomes smaller compared that rigid body. 

Since was not possible correlate the results the experiments 
the motion bubbles the gravity pressure field terms non-dimensional 
parameters formed from the usual liquid properties (viscosity, surface ten- 
sion, density), further work freely rising bubbles necessary before the 
results obtained from such tests can utilized and before the more com- 
plicated behavior bubbles variable pressure gradients can under- 
stood. Particularly, understanding the reason for the transition bub- 
bles from fluid ‘rigid’ bodies well criterion for this transition 
most desirable, since such transition might influenced the magnitude 
the pressure gradient. the region Reynolds numbers where the bub- 
bles behave like rigid spheres, the pressure gradient probably has effect 
the drag coefficient. 


SUMMARY 


the size the bubbles was increased the tests, change bubble 
shape from spherical ellipsoidal spherical cap shape was observed 
all liquids. The volumes which these transitions occur, however, varied 
with the properties the liquid. For spherical bubbles given volume the 
results show that the viscosity the liquid the most important property 
determining the rate rise. For ellipsoidal bubbles the surface tension 
assumes greater importance. Spherical cap bubbles rise independently 
the properties the liquid. 

The results show that the motion air bubbles rising their terminal 
velocity gravity field cannot described completely use dimen- 
sionless parameters formed from the usual liquid properties (viscosity, sur- 
face tension, density), the equivalent radius the bubble, the acceleration 
gravity, and the terminal velocity. 

The drag coefficients tiny spherical bubbles coincide with those cor- 
responding rigid spheres. With increase bubble size, decrease the 
drag compared that rigid spheres occurs some liquids. This 
change the drag due the development circulation inside the bubble. 
The drag curves the spherical bubbles rising various liquids fall be- 
tween two limiting curves, namely the drag curves rigid and fluid spheres, 
respectively. was not possible determine criterion for the transition 
region the bubbles from ‘rigid’ fluid spheres. 

Ellipsoidal bubbles occur for different ranges Reynolds numbers for 
the various liquids. For liquids low number (say less than 
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minimum the drag curve reached Reynolds numbers the order 

250. These minima occur near the transition from spherical ellipsoidal 
shape. Such minima the drag curve are not obtained for liquids high 

number. For the liquids used, transition spherical caps completed 

Weber number about twenty. 

The drag spherical cap bubbles are independent bubble 
size and have constant value 2.6. The rate rise these bubbles 
function the equivalent radius given the experimentally determined 
relation: 


For bubbles (ranging equivalent radius from 0.03 0.25 cm) rising 
tap water, increased drag compared bubbles clean (filtered 
distilled) water was observed. The presence certain surface-active sub- 
stances the water similarly increases the drag bubbles (ranging 
equivalent radius from 0.03 0.30 cm) compared bubbles pure wa- 
ter. Beyond certain critical concentration these surface-active sub- 
stances, increase concentration has relatively little influence the 
drag the bubbles. 

Tests determine the effect the container walls the velocity rise 
indicate the absence such effect for the range bubble volumes and con- 
tainer sizes tested. 
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APPENDIX LIST SYMBOLS 


Drag coefficient (8/3)g 


Drag bubble 

Coefficient surface viscosity 
Acceleration due gravity 
Length parameter bubble 


Dimensionless parameter 


Equivalent radius bubble, i.e. radius 
sphere equal volume 


Reynolds number 


Reynolds number inside bubble 
Terminal velocity 


Coefficient dynamic viscosity liquid medium 
Coefficient dynamic viscosity gas 

Density liquid medium 

Density gas 


Surface tension 


(8) 
(9) 
(10) 
(11) 
(12) 
(13) 
(14) 
p! 
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TABLE 


Summary Liquid Properties 


Surface 

Temperature Viscosity Density Tension 

deg poises dynes/cm 


Water 72.9 0.26 
Water 72.6 0.24 
Cold Water 74.8 1.08 
Hot Water 68.1 0.307 
Solution 32.8 2.78 
Oil 20.7 1.45 
Methy! Alcohol 21.8 0.89 


percent Corn Syrup 79.2 0.155 
and Water 


percent Corn Syrup 79.9 
anu Water 


percent Glycerine 
and 


percent Glycerine 
and Water 


percent Ethy! Alcohol 10% 
Water 


Olive Oil 0.716 
Syrup 0.92 10° 


TABLE 
Bubble Shapes Function Reynolds Number 
Reynolds Number 


than 400 5000 larger than 5000 
Cold water 275 3000 


1500 


Methyl Alcohol 4000 


Corn Syrup 
and Water 
68% Corn Syrup 
and Water 


| | | 
| | | 
4000 
| 
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Equivalent Radius in centimeters 


Comparison Terminal Velocities Air Bubbles Various 
Liquids 
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Fig. Typical Paths Bubbles Rising Various Liquids 
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